The origin of the stabilized simple-cubic structure in polonium 
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The origin of the stabilized simple-cubic (SC) structure in Po is explored by using the first 
principle band calculations. We have found that the prime origin is the inherent strong spin-orbit 
(SO) interaction in Po, which suppresses the Peierls-like structural instability as usually occurs in 
p-bonded systems. Based on the systematic analysis of electronic structures, charge densities, Fermi 
surfaces, and susceptibilities of Se, Te, and Po, we have proved that the stable crystal structure 
in VIA elements is determined by the competition between the SO splitting and the crystal field 
splitting induced by the low-symmetry structural transition. The trigonal structure is stabilized in 
Se and Te by the larger crystal field splitting than the SO splitting, whereas in Po the SC structure 
is stabilized by the large SO splitting. 

PACS numbers: 61.66.Bi, 71.70.Ej, 71.20.Gj 



Among the elements in the periodic table, only Po 
(polonium) is known to crystallize in the simple-cubic 
(SC) structure in nature 1 . Radioactive Po, which was 
discovered by Pierre and Marie Curie in 1898, requires 
difficult sample preparation, and so not many experi- 
mental and theoretical studies have been carried out. Po, 
which is located in the VIA column of the periodic table, 
has a valence electron configuration of 6s 2 6p 4 (2=84), 
and exhibits two metallic phases of a and j3 which are 
formed in the SC and the rhombohedral (trigonal) struc- 
ture below and above ~ 348-ftT, respectively 2 ^. Note that 
the isoelectronic elements located in the same VIA col- 
umn, Se (selenium) and Te (tellurium) , crystallize in the 
trigonal structure (space group P3i21) with the helical 
chain arrangements of atoms, which run parallel to the 
crystallographic c axis of the hexagonal setting. Hence 
the coordination number is two in trigonal Se and Te in 
contrast to six in simple-cubic Po (SC-Po). The ratio 
of the intrachain (di ) bond length to the interchain (^2) 
bond length is & = 1.59, 1.23, and 1.0 for Se, Te, and Po, 
respectively. The coordination number of two in trigonal 
Se and Te can be understood based on the simple octet 
rule^. 

The reason why metallic Po has a stable SC struc- 
ture has not been fully addressed yet. By using the rel- 
ativistic parameterized extended Hiickel method, Lohr 5 
has found a hint that the helical structure as a distor- 
tion of a SC structure might be quenched in the case 
of Po due to its very large spin-orbit (SO) coupling. 
However, their calculation was not self-consistent and 
adopted the atomic value of the SO coupling parame- 
ter. They did not study the structural energetics. Re- 
cently, using the ab-initio pseudopotential band method, 
Kraig et alA showed that the SC structure is preferred 
by Po to face-centered or body-centered cubic structure. 
They argued that the large s-p splitting in Po would pro- 
duce the stable SC structure. However, they could not 
explain why this happens only in Po but not in other el- 
ements with similarly large s-p splittings. Moreover, the 
SO interaction was not taken into account in their band 



calculations. Lach-hab et al£ studied the structural ener- 
getics for Po using the tight-binding (TB) band method. 
After having determined the TB parameters by fitting 
the TB bands to the semirelativistic full-potential lin- 
earized augmented plane- wave (FLAPW) 9 bands, they 
took account of the SO effect a posteriori by employing 
the atomic value of the SO coupling parameter. They 
found that the SC structure is the most stable among the 
close-packed structures for both cases with and without 
the SO effect included. They did not consider the lower 
symmetry structures with the internal atomic relaxation, 
as observed in Se and Te. 

A systematic crystallographic transformation can oc- 
cur from the SC lattice to the trigonal lattice. As shown 
in Fig.QJa), the SC lattice can be described by a hexag- 
onal lattice in which the hexagonal a axis corresponds to 
the face diagonal and the hexagonal c axis to the body 
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FIG. 1: (a) A hexagonal unit cell with tripled volume of a SC 
unit cell, (b) The hexagonal Brillouin zone (BZ) compared to 
the SC BZ. Note that T-A and T-L in the hexagonal BZ cor- 
respond to ^P-R and T-X in the SC BZ. (c) Schematic total 
energy curves for Se, Te, and Po in the hypothetical struc- 
tural space. T and SC represent trigonal and simple-cubic 
structures, respectively. Both in the GGA and GGA+SO, 
trigonal Se and Te are more stable. For Po, trigonal Po is 
more stable in the GGA, while SC-Po becomes more stable 
in the GGA+SO. 
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FIG. 2: (a)-(b) The GGA and GGA+SO band structures of 
SC-Po. (c)-(d) DOSs of SC-Po and trigonal Po. The latter 
corresponds to the stable phase obtained after the geometry 
optimization in the GGA scheme. In each frame, DOSs ob- 
tained by the GGA and GGA+SO schemes are compared. 
(e)-(f) The GGA and GGA+SO band structures of SC-Po 
drawn in the frame of the hexagonal BZ. There are several 
band crossings at Ef in the GGA, but not in the GGA+SO. 



diagonal of cube. Then the size of the unit cell becomes 

tripled with c/a — y|. The trigonal structure of /3-Po 

corresponds to the elongated SC structure simply along 
the body diagonal direction. In the trigonal structures of 
Se and Te, the internal atomic position parameters are 
also changed in addition to elongation^. 

To explore the origin of stabilized SC structure in Po, 
we have investigated the electronic and structural prop- 
erties of Se, Te, and Po systematically by using the first- 
principles FLAPW band method incorporating the SO 
interaction. For the calculations, the QMD-FLAPW 9 
and the WIEN2K 10 codes are utilized. The band struc- 
ture results obtained from both codes are qualitatively 
coincident with each other. Below we will use mainly 
WIEN2K results, since the local orbital basis is im- 



plemented in WIEN2K in order to describe the rela- 
tivistic 6pi/2 wave function correctljsii. The exchange- 
correlation interaction is treated by the generalized gra- 
dient approximation (GGA) 12 - 13 . 

We have first studied the structural energetics by op- 
timizing the atomic positions for both trigonal and SC 
structures of Se, Te, and Po with (GGA+SO) and with- 
out (GGA) the SO interaction included. In the GGA 
scheme, the trigonal structures are more stable than SC 
structures by 0.092, 0.031, and 0.004 eV/atom for Se, 
Te, and Po, respectively. The GGA+SO scheme, how- 
ever, yields different results. The incorporation of the SO 
coupling reduces the stability of the trigonal structure. 
For Se and Te, the trigonal structure is still more sta- 
ble than the SC structure by 0.090 and 0.021 eV/atom, 
respectively. In contrast, for Po, the SC structure be- 
comes more stable than the trigonal structure by 0.016 
eV/atom. Thus the total energy behaves as in Fig.^c). 
There are two local minima in the hypothetical struc- 
tural space, the trigonal (T) and the SC structure. The 
SO coupling effect does not change the global minima in 
Se and Te, but changes that in Po from trigonal Po to 
SC-Po. This finding indicates that the SO interaction 
plays an essential role in stabilizing the SC structure of 
Po. 

Indeed band structures in Fig.|2Ia)-(b) reveal that the 
effect of the SO interaction is substantial for SC-Po, as 
compared to that for Se and Te. The lowest bands in 
Fig-EIa)-(b) correspond to Po 6s states, while the bands 
near the Fermi level {Ep) to Po 6p states. Due to the 
large energy gap between 6s and 6p states, there will 
not be a considerable mixing between them, and thus 
Po can be classified as a p-bonded system. When the SO 
interaction is taken into account, the degenerate 6p bands 
along T-X, T-M, and T-R are split so that the doubly 
degenerate 6pi/2 band is separated from the upper half- 
filled 6p 3 / 2 band. Due to the SO splitting, the small hole 
Fermi surface at T disappears. The SO splitting of p 
states at T for SC-Po is as much as 3.16 eV, which is 
much larger than those for SC structures of Se and Te, 
0.53 and 1.05 eV. 

In Fig.|2Ic)-(d), we have compared the density of states 
(DOS) of SC-Po and that of trigonal Po, which are ob- 
tained by the GGA and GGA+SO schemes. Trigonal 
Po here corresponds to the stable phase determined by 
the geometry optimization in the GGA scheme. Notable 
in the GGA-DOS is the pseudogap feature at Ep for 
trigonal Po, which would give rise to the semi-metallic 
behavior. The pseudogap appears between the nonbond- 
ing lone-pair and the antibonding Po 6p state, because of 
the splitting between the bonding and antibonding Po 6p 
state in the chain arrangement of trigonal Poii Hence 
the DOS at Ep, N(Ep), becomes much reduced in trig- 
onal Po from that in SC-Po. In fact, the stable trigonal 
phases of Se and Te result from this gap formation near 
Ep. In the case of Se, the energy splitting between the 
bonding and the antibonding state is so large that the 
system becomes insulating, while Te is semimetallic. 
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FIG. 3: The occupied charge density distributions of Po 6p 
bonding states in the energy range of —6.0 ~ —3.0 eV. (a) SC- 
Po (GGA), (b) trigonal Po (GGA), (c) SC-Po (GGA+SO), 
and (d) trigonal Po (GGA+SO). 



In the GGA+SO scheme, the &P1/2 state is separated 
out, and so the band shape and N(Ep) are changed a 
lot. Occupied band widths are broadened with more 
weight at low energy side, which produces energy gains 
for both SC-Po and trigonal Po. The stable SC-Po in 
the GGA+SO scheme reflects that the SO induced en- 
ergy gain is larger for SC-Po than for trigonal Po. It is 
because N(Ep) for trigonal Po is rather enhanced pro- 
ducing some energy loss. Note that, for Se and Te, both 
the GGA and GGA+SO schemes yield similar DOS, im- 
plying no noticeable energy gain from the SO interaction. 
These results suggest that the stable crystal structure in 
VIA elements is mainly determined by the competition 
between the SO splitting and the crystal field splitting 
due to the low symmetry structure. That is, for Se and 
Te, the larger crystal field splitting between the intra- 
chain bonding and antibonding states stabilizes the trig- 
onal structure, whereas, for Po, the larger SO splitting 
stabilizes the SC structure. 

The stable trigonal structure of Se and Te can be un- 
derstood in terms of the three dimensional (3D) Peierls 
distortion of the SC structur e) 4 ! 8 ; 15 . In general, p-bonded 
systems would favor the SC structure to maximize the 
directional p a bonding. However, three mutually orthog- 
onal linear chains in the SC structure will be easily dis- 
torted by the Peierls mechanism when the linear chain is 
partially filled. In the case of Se and Te, each chain is 
2/3-filled, and so the energy gain can be achieved by the 
trimcrization with the short-long-long bond alternation. 
Decker et ali^> recently demonstrated this mechanism for 
Te by comparing the band structures of the undistorted 
SC and the distorted trigonal structure. The structural 
distortion from the SC to the trigonal structure induces 
a splitting of degenerate bands near Ep, which yields the 
energy gain for trigonal Te. 

We have examined the above scenario for SC-Po. Fig- 
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FIG. 4: Fermi surfaces of SC-Po derived from the third 
(green) and fourth (blue) bands in the GGA (a) and the 
GGA+SO (b). (c) Susceptibilities of SC-Po in the GGA and 
the GGA+SO. 



ureEIe)-(f) display the GGA and GGA+SO band struc- 
tures of undistorted SC-Po which are projected in the 
hexagonal BZ. The GGA band of Fig.[3fe) shows several 
band crossings at £f which are created by the back- 
folding of bands of the larger SC BZ into the smaller 
hexagonal BZ (see Fig.^b)). Then, as in Tei£, a Peierls- 
like structural distortion would split the band crossings at 
Ep to yield the energy gain. In contrast, the GGA+SO 
band of Fig.^f) shows the band crossings not at Ep but 
at higher or lower energy side far apart from Ep. This 
is because of the substantial band shift by the SO split- 
ting. Accordingly, no energy gain would result from the 
Peierls-like structural distortion. This proves how the 
SO interaction stabilizes the SC structure of Po. When 
Decker et alJ£° demonstrated this mechanism for Te, they 
did not take into account the SO interaction. The SO in- 
duced band shift in Te, however, is not large enough to 
suppress the Peierls-like instability. 

The bonding characters can be investigated by the 
charge density plot. Figure shows the occupied charge 
density distributions of Po 6p bonding states in the en- 
ergy range of —6.0 ~ —3.0 eV for both SC-Po and trigo- 
nal Po. They are plotted on the special plane of hexago- 
nal unit cell to see the bonding character along the helical 
zig-zag chains in the c-direction. The GGA charge den- 
sity for SC-Po (Fig. Ola)) shows larger bonding charac- 
ter between Po atoms than the GGA+SO charge density 
(Fig.[3fc)). This implies that the directional bonding be- 
tween atoms becomes weakened when incorporating the 
SO interaction, and so the charge density becomes more 
isotropic. This feature is more clearly seen in trigonal 
Po which manifests the prominent chain nature. The 
GGA charge density for trigonal Po (Fig.[3fb)) shows a 
pronounced intrachain bonding character along the zig- 
zag chains. In the GGA+SO charge density (Fig. Eld)), 
the SO interaction weakens the intrachain directional 
bondings so that the anisotropic chain nature is sup- 
pressed. This is consistent with the enhanced N(Ep) 
in the GGA+SO DOS in Fig.^d). The charge density 
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plots reveal that the SO interaction weakens the direc- 
tional bondings between Po atoms and so suppresses the 
ID chain nature and the corresponding Peierls-like struc- 
tural instability. 

The structural transition can also be studied by an- 
alyzing the behavior of the charge susceptibility Xo(o)- 
Xo(q)' s m Fig- EI C ) are obtained from the band struc- 
tures of SC-Po. Xo(q) in the GGA scheme shows the 
highest peak near Q = §r-R, indicating a possible struc- 
tural instability at this Q. This Q vector is coincident 
with the reciprocal lattice vector connecting the hexago- 
nal BZ boundaries including A (see Fig. db)), at which 
the bands of SC-Po along T-R are folded back. There- 
fore, this peak is closely related to the trigonal distortion 
of the SC structure along the (111) direction. Notewor- 
thy in Fig. 2Jc) is the substantial reduction of this peak 
°f Xo(q) when incorporating the SO interaction. The ab- 
sence of the structural instability in Po would be closely 
correlated to this behavior of Xo(<z)- One can also note 
the peaks in xo(q) near q = |r-X, fr-M, |X-M. These 
peaks and the peak near Q are expected to arise from the 
quasi-lD Fermi surfaces of SC-Po, which are formed by 
the third and fourth bands in Fig. Ufa). Namely, each q 
producing the peak in xo(o) corresponds to the nesting 
vector of the quasi-lD Fermi surfaces projected in each 
symmetry plane. The GGA Fermi surfaces in Fig. Ufa), 
which are drawn in the (220) plane of the SC BZ, exhibit 
the clear quasi-lD nature with the corresponding nesting 
vectors. As shown in Fig.HJb), the SO interaction breaks 
these quasi- ID Fermi surfaces into pieces, and so reduces 
the nesting effect. This situation in Po is different from 
that in Se and Te, in which the ID nature is preserved 



even with the SO interaction, so that the nesting effect 
is still active. 

Finally, it is worthwhile to discuss the stable structures 
of neighboring elements of Po, Bi (Z=83) and At (Z=85), 
both of which are also p-bonded metallic systems with 
the large SO interaction. In nature, Bi crystallizes in 
the trigonal (a-arsenic) structure^. As for radioactive 
At, no structural information is available because of its 
too short half-life of only 8.3 hours, a-arsenic structure 
of Bi can be understood as arising from the Peierls-likc 
distortion, as discussed abov o 4 i 15 . Since the 6p-bands in 
Bi are half-filled, the tendency of the Peierls instability 
is predominating so as to drive the structural distortion 
despite the large SO interaction. In the case of At, the 
6p- bands are 5/6-filled so that the SO interaction would 
dominate over the Peierls instability as in Po. Hence the 
SC structure is expected to be stabilized for At. 

In conclusion, we have clarified that the origin of the 
stabilized SC structure in Po is its inherent strong SO 
interaction. The large SO interaction in Po weakens the 
directional bondings between atoms so as to suppress the 
Peierls-like distortion. Our study also reveals that the 
stable crystal structure in VIA elements is determined by 
the competition between the SO splitting and the crystal 
field splitting induced by the structural transition to a 
lower symmetry. For Se and Te, the larger crystal field 
splitting between the intra-chain bonding and the anti- 
bonding state stabilizes the trigonal structure, whereas, 
for Po, the larger SO splitting stabilizes the SC structure. 
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